Self-consistency of relativistic observables with general relativity in 
the white dwarf-neutron star binary pulsar PSR J1141-6545. 
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ABSTRACT 

Here we report timing measurements of the relativistic binary pulsar 
PSR J1141— 6545 that constrain the component masses and demonstrate that 
the orbital period derivative Py, = (—4 ± 1) x 10~^^ is consistent with gravita- 
tional wave emission as described by the general theory of relativity. The mass 
of the neutron star and its companion are 1.30±0.02 Mq and 0.986±0.02 
respectively, suggesting a white dwarf companion, and extending the range of 
systems for which general relativity provides a correct description. On evolu- 
tionary grounds, the progenitor mass of PSR J1141— 6545 should be near the 
minimum for neutron star production. Its mass is two standard deviations below 
the mean of the other neutron stars, suggesting a relationship between progenitor 
and remnant masses. 

Subject headings: gravitational radiation, binary pulsars, neutron star masses, 
general relativity, PSR J1141-6545 

1. Introduction 



Our knowledge of neutron star masses (Thorsett and Chakrabarty 1999) is derived largely 
from three relativistic binary pulsars with eccentric orbits and probable neutron star compan- 
ions (Weisberg and Taylor 2003; Stairs et al. 2002; Deich and Kulkarni 1996). The masses 
of all known neutron stars exhibit a remarkably narrow scatter (O.O4M0) about their mean 
of 1.35 Mq. If one restricts attention to the neutron star pairs, the mean is 1.37 Mq. These 
binary pulsars have also been pivotal in tests of the general theory of relativity (Taylor, 
Fowler and McCulloch 1979). 

Gravitational waves have never been observed directly. Our best quantitative evidence 
for their existence comes from observations of the binary pulsar PSR B1913+16 (Damour and 
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Taylor 1991). This pulsar has a negative orbital period derivative that is consistent with the 
emission of gravity waves at the 0.4% level. Unfortunately, the binary PSR B1534+12 is so 
close to the Earth that its proper motion induces a time- varying Doppler term into the timing 
which prevents separation of the relativistic contribution to the orbital period derivative from 
the kinematic one (Bell and Bailes 1996; Stairs et al. 2002). Until an independent distance 
can be established, this pulsar's orbital period derivative cannot verify the predictions of the 
theory to better than ~ 20%. Nevertheless, measurement of PSR Blo34+12's post-Keplerian 
parameters have provided other tests of general relativity (Stairs et al. 2002), and binary 
pulsars such as PSR J0437-4715 have proved that the theory is consistently in agreement 
with the most precise observations (van Straten et al. 2001). 

Discovered in the Parkes multibeam survey PSR J1141— 6545 is a 4.7*^ binary pulsar 
with an eccentricity of 0.17 that consists of a 394 ms relatively youthful pulsar of characteris- 
tic age ~1.4 Myr and an unidentified companion star. Early observations (Kaspi et al. 2000) 
of the periastron advance suggested a pulsar mass Mp < 1.35 Mq, and a companion 
Mc > 0.97 M0. The observed tight distribution of neutron star masses near 1.35±0.04 
Mq where the quoted error refiects a Gaussian fit to the observed distribution, suggested 
that this system probably consisted of a Mp ~ 1.3 M0 neutron star and Mc ~ 1 Mq white 
dwarf companion. PSR J1141— 6545 was therefore unique among the relativistic binary 
pulsars. Its component masses argued for a white dwarf companion yet its orbit was highly 
eccentric making it the first of a new class of object. To date no white dwarf has been 
optically identified. 

PSR ,11141—6545 is an exciting target for timing. Firstly, scintillation observations 
have demonstrated a clear modulation of the timescale of the interstellar "twinkling" allow- 
ing derivation of the system's angle of inclination and runaway velocity (Ord, Bailes and van 
Straten 2002a) . Two degenerate fits to the changing scintillation timescale were possible. 
The observed distribution of neutron star masses favoured the solution with an inclination an- 
gle (i = 76±2.5°) and center of mass velocity in the plane of the sky of y ~ 115 km s~^. These 
results provide an independent estimate of the system orientation to be compared with the 
timing. Secondly, neutral hydrogen observations (Ord, Bailes and van Straten 2002b) have 
demonstrated that this pulsar is in a favourable location in the Galaxy where the observed 
orbital period derivative will, as in the case of PSR B1913 16, be dominated by the general 
relativistic contribution. Finally, contrary to PSR B1913 I 16, which is a very symmetric sys- 
tem with two neutron stars of nearly equal mass, PSR J1141— 6545 is a very dissymmetric 
system as it comprises a strongly self-gravitating neutron star and a (relatively) weakly self- 
gravitating white dwarf. Timing this system therefore tests new aspects of relativistic theo- 
ries of gravity (Damour and Esposito-Farese 1992; Taylor et al. 1992; Arzoumanian 2003). 
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2. Observations and Data Analysis 



Observations of PSR J1141— 6545 were conducted from January 2001 until February 
2003 at the Parkes 64 m radio telescope as part of a scintillation study. The data from two 
orthogonal linear polarisations were detected and summed using a multi-channel filterbank 
consisting of 2x512x0.5 MHz filters. These data were 1-bit sampled every 125 fis and 
written to magnetic tape for offline processing on the Swinburne University of Technology 
supercomputer. Observations were usually greater than one orbit in duration, and there 
were ten sessions, nine of which took place in the 15 months prior to Feb 2003. The central 
observing frequency was 1390 MHz and total intensity pulse profiles were formed by folding 
the raw data at the topocentric period of the pulsar for integration times of 60 s. 

A fit to the arrival times was made for the standard pulsar spin, astrometric, and Kep- 
lerian terms as well as three post-Keplerian parameters using a relativistic model (Damour 
and Deruelle 1986). After incorporation of archival timing data from the Australia Tele- 
scope National Facility archives from July-August 1999 we found that our measurements of 
the periastron advance uj, transverse Doppler and gravitational redshift, 7 term, and orbital 
period derivative Pb became much more significant. These values are displayed in Table 
1. The archival data consisted mainly of data taken with a 2x96x3 MHz filterbank sys- 
tem. We fit for a single arbitrary offset between the two datasets to allow for the change in 
instrumentation. 



Newton's laws provide a relation between the system's masses and inclination angle i, 
and the common observables orbital period Pb and projected semi-major axis x = {asmi)/c. 
This is often referred to as the pulsar mass function /p, 



here M is the sum of the masses, G is Newton's gravitational constant and c is the speed of 
light. Mp and Mc refer to the mass of the pulsar and companion star respectively. 

General Relativity allows us to use the post-Keplerian terms to derive component 
masses. The advance of periastron gives us the sum of the masses. 



3. Interpretation 




(1) 




(2) 



Here e is the orbital eccentricity. The combined transverse Doppler and gravitational redshift 
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term gives a different relation between tlie masses. 



(3) 



General relativity predicts that gravitational wave emission from binary systems removes 
energy that results in an orbital decay. The system is over-determined if one can also 



Thus, binary pulsars such as PSR J1141— 6545 can be used to test the validity of the general 
theory of relativity. The allowed range of masses for the system are shown in Figure 1. 
Other post-Keplerian parameters such as the range and shape of the Shapiro delay were not 
significant in the fit. 

The companion is almost certainly a white dwarf, unless neutron star masses extend 
over a very large range. The mean and standard deviation of the 7 neutron stars shown in 
Figure 2 would place the companion of PSR J1141— 6545 seven standard deviations from the 
mean of the others which appears unlikely. 

Young pulsars often exhibit significant timing noise which make timing proper motions 
unreliable on the timescale of a few years. A formal fit to the limited data yielded proper 
motions in both right acension and declination of -12±3 milli-arc seconds per year. Whether 
the apparently 4-sigma proper motion is real or not will become apparent with continued 
timing. For now we are treating the result with some caution. Timing noise has a much 
longer period (months to years) than the orbital period of our pulsar (4.7h), and so although 
the timing proper motion is potentially non-physical, the derived orbital parameters are not 
subject to contamination by timing noise. 

The self-consistency of the system with general relativity is very pleasing. The orbital 
period derivative is now a 4-sigma result. The relative error in this term scales as the 
length of the observing span so by the end of the decade it should be near 1%. 

PSR J1141— 6545is therefore likely to be a very important new testing ground for general 
relativity. The dissymetry of the system tests different aspects of the theory since one of the 
stars is not strongly self-gravitating. Its favourable location in the Galaxy means that the 
measured orbital period derivative is close to the intrinsic value (Ord, Bailes and van Straten 
2002a). Unfortunately, the pulsar is young, and unlikely to be as intrinsically stable as PSR 
B1913-rl6, which has a much weaker magnetic field strength. Whether this is important or 
not for post-Keplerian orbital parameters remains to be seen. In any case, it will be difficult 



measure the orbital period derivative Pb- 





(4) 
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to surpass PSR B1913+16's relative errors on terms such as the relativistic contribution to 
the orbital period derivative. 

4. Discussion 

Large-scale simulations of binary pulsar progenitors predict a large population of ec- 
centric pulsar- white dwarf binaries (Tauris and Sennels 2000), albeit with a short observ- 
able lifetime. PSR J1141— 6545is thought to arise from the evolution of a binary in which 
both stars are initially less than the critical mass Merit required to produce a supernova 
(Tauris and Sennels 2000). The initially more-massive star transfers mass to its companion 
before becoming a white dwarf. If sufficient matter can be accreted by the initially lower- 
mass star it will exceed Merit and produce a pulsar. Should the system remain bound, we 
will be left with a young neutron star orbiting a white dwarf companion in an eccentric 
orbit. PSR J1141— 6545's mass Mpis therefore very important, because we know that the 
progenitor must have been less than 2 Merit - Mc . This is the first time we have an upper 
limit to the progenitor mass of a neutron star, albeit in terms of Merit • All previous neutron 
star mass measurements are from systems where both components were large enough to 
form neutron stars, which makes the progenitor masses far more difficult to constrain as our 
theoretical understanding of the mass of neutron star progenitors is limited. 

Radio pulsar masses cannot be easily determined. Neutron stars can theoretically exist 
over a range of masses (Timmes, Woosley and Weaver 1996) that is much broader than that 
inferred from observations of radio pulsars ( Weisberg and Taylor 2003, Stairs et al.2002, 
Deich and Kulkarni 1996). Until these measurements, the mean and standard deviation of 
the six accurately known neutron stars was just 1.37±0.04 Mq. PSR J1141— 6545's mass 
of 1.30±0.02 is the smallest yet determined with any accuracy. It's mass is two standard 
deviations below the mean of the others as shown in Figure 2. 

In the absence of any asymmetry in the explosion, component masses and orbital parameters 
define the progenitor system mass and runaway velocity (Radhakrishnan and Shukre 1986). 
If the explosion were symmetric, the pre-supernova star in the PSR J1141— 6545 system 
would have been only Mpre = M(l + e) — Me = 1.7 M0. It seems unlikely that the pre- 
supernova star could be such a small fraction of its original mass > Merit ~ 8 — 11 Mq. 
Thus PSR J1141— 6545's orbital configuration is therefore a strong argument in favour of 
asymmetric supernova kicks (Dewey and Cordes 1987; Bailes 1989). Such kicks can misalign 
the spin and orbital angular momentum vectors of the pulsar, allowing detectable precession 
due to spin orbit coupling(Damour and Ruffini 1974; Barker and O'Connell 1975). This has 
been observed in the binary pulsar PSR 1913+16 (Kramer 1998). In PSR J1141-6545 the 
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timescale for complete precession is only ~ 200 yr. Its high orbital inclination could lead to 
significant changes in our line of sight to the pulsar's spin axis on the timescale of a decade 
or so. This might explain the non-detection of the pulsar in previous surveys of the region 
and we would expect rapid evolution of the pulse profile as suggested by (Kaspi et al. 2000) . 

PSR J1141— 6545 is a very significant object. It has already validated scintillation 
speed methodologies, provided the first estimate of a neutron star mass where we have 
firm constraints on the progenitor, proved that mass transfer can transform a star that 
would otherwise become a white dwarf into a neutron star and extended the range of 
systems for which general relativity provides an accurate description. Its orbital param- 
eters suggest that the explosion was probably asymmetric. As the errors in orbital pe- 
riod derivative decrease, this system will test new aspects of theories of relativistic grav- 
ity that were impossible with the other relativistic binaries due to its unique composition 
(Taylor et al 1992; Damour and Taylor 1992). Together, PSR B1913+16, PSR B1534+12 
and PSR Jl 141-6545 will strongly constrain relativistic theories of gravity. 
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Table 1 Observed and derived parameters 



Right ascension, a (J2000) 


llMl'"0^022(6) 


Declination, 5 (J2000) 


-65°45'19'.'089(9) 


Pulse period, P (ms) 


393.8978340370(2) 


Reference epoch (MJD) 


51369.8525 


Dispersion Measure cm^ pc~^ . . 


116.048(2) 


Period derivative, P (10^^) . . . 


4.294593(3) 


Orbital period, P^ (days) 


0.1976509587(3) 


Projected semi-major axis x (s) 


1.85894(1) 


Orbital eccentricity, e 


0.171876(2) 


Epoch of periastron. To (MJD) . 


51369.854553(1) 


Longitude of periastron, cu (°) •• 


42.457(2) 


Number of arrival times 


7180 


RMS Residual (/xs) 


114 


Damour-Deruelle PK parameters 




7 (s) 


0.00072(3) 


u (°yr-i) 


5.3084(9) 


Pb(10-^^) 


-0.43(10) 


Damour-Deruelle GR fit 




Companion mass, Mc (Mq) 


0.986(20) 


Pulsar mass, Mp (M©) 


1.30(2) 


Sum of masses, M (Mq) 


2.2883(5) 


Derived Parameters 




Orbital inclination, i (°) 


i > 75 



Best-fit physical parameters and their formal la er- 
rors were derived from arrival time data by min- 
imizing an objective function, as implemented in 
TEMPO (http://pulsar.princeton.edu/tempo). The 
timing model used is TEMPO'S DD relativistic bi- 
nary model (Damour and Deruelle 1986) which al- 
lows us to fit for Lu, Pb and 7 separately. To de- 
rive the best estimate of the companion and pulsar 
masses, a fit to only M and was performed using 
tempo's DDGR model. This model assumes gen- 
eral relativity is correct. Parenthesized numbers rep- 
resent the l-cr uncertainty in the last digits quoted, 
and epochs are specified using the Modified Julian 
Day (MJD). 
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Fig. 1. — Constraints on the pulsar and companion masses. The shaded area is excluded 
from a consideration of Kepler's laws. The sum of the masses M = 2.2883 ± O.OOO6M0 
is very accurately defined from the relativistic advance of periastron, uj, and is shown by 
the straight solid line with a slope of -1. The relativistic term 7 constrains the system to 
lie between the dashed lines, and the measured orbital period derivative A suggests that 
the system masses lie between the pair of curved solid lines. This implies an inclination 
angle i > 75°, consistent with constraints from scintillation studies i = 76 ± 2.5°. The 
derived component masses suggest that the system consists of a 0.986±0.020 Mq white 
dwarf and a 1.30±0.02 Mopulsar. If the system radiates energy at the rate predicted by 
general relativity, the measured orbital period derivative Pb = (— 4±1) x 10"^'^ is consistent 
with the predicted value of (—3.8 x 10^^) derived from the component masses. 
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Fig. 2. — Neutron star masses derived from relativistic binary pulsars. Error bars are shown 
for all data, but in the case of PSR B1913 + 16 and PSR B1534+12 they are too small to be 
visible. Masses are taken from Deich and Kulkarni (1996), Stairs et at 2002 and Weisberg 
and Taylor (2003) . PSR J1141 —6545 is the first neutron star not in a neutron star-neutron 
star binary to have its mass so precisely determined. On evolutionary grounds we know 
its progenitor was near the lower end of the mass distribution for pulsar progenitors. The 
low remnant mass suggests that the neutron star mass is related to that of the progenitor. 



